Separation and identification of dominant mechanisms in double photoionization 
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Double photoionization by a single photon is often discussed in terms of two contributing mech- 
anisms, knock-out (two-step-one) and shake-off with the latter being a pure quantum effect. It is 
shown that a quasi-classical description of knock-out and a simple quantum calculation of shake-off 
provides a clear separation of the mechanisms and facilitates their calculation considerably. The 
relevance of each mechanism at different photon energies is quantified for helium. Photoionization 
ratios, integral and singly differential cross sections obtained by us are in excellent agreement with 
benchmark experimental data and recent theoretical results. 
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PACS numbers: 3.65.Sq, 32.80.Fb, 34.80.Dp 

Our understanding of dynamical processes often rests 
on isolating approximate mechanisms which leave charac- 
teristic traces in the measured or computed observables. 
A prime example is double photoionization. After the 
initial absorption of the photon by the primary electron 
the subsequent redistribution of the energy among the 
electrons is often discussed in terms of two mechanisms 
0, knock-out (KO) (sometimes called 'two-step-one' ||) 
and shake-off (SO) || f|, ||. The first mechanism de- 
scribes the correlated dynamics of the two electrons as 
they leave the nucleus where the primary electron has 
knocked out the secondary electron in an (e, 2e)-like pro- 
cess. The second mechanism accounts for the fact that 
absorption of the photon may lead to a sudden removal 
of the primary electron. This causes a change in the 
atomic field so that the secondary electron relaxes with a 
certain probability to an unbound state of the remaining 
He + ion, i.e., the secondary electron is shaken off. 

Apart from general properties, e.g., the prevalence of 
shake-off at high photon energies, it is difficult to separate 
the processes. However, they are distinct with respect 
to their quantum nature: shake-off is a purely quantum 
mechanical phenomenon while knock-out dynamics oc- 
curs classically as well as quantum mechanically. This 
opens up the possibility to separate shake-off and knock- 
out by calculating the latter (quasi-)classically, provided 
the quasi-classical approximation to knock-out is good. 
Clearly, the quasi-classical KO mechanism does not con- 
tain any part of SO (which is purely quantum). 

The two phases of double photoionization, initial ab- 
sorption and redistribution of the energy, can be ex- 
pressed by the relation 



<4 + = <W£ + 



(1) 



where X stands for either shake-off or knock-out. In 
the following, we evaluate Px + f° r full fragmentation 
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of the ground state and use the experimental data of 
Samson et al. || for a a ^ s . We obtain the classical dou- 
ble escape probability for KO with a classical-trajectory 
Monte-Carlo (CTMC) phase space method. CTMC has 
been frequently used for particle impact induced frag- 
mentation [0, H |l0| with implementations differing 
typically in the way the phase space distribution, p(T), 
of the initial state is constructed. Details of our approach 
will be published elsewhere, here we summarize the im- 
portant steps only. 

Within our phase space approach the double escape 
probability ijfo ^ s formally given by 



P 



KO 



lim / dr 



- -p++ 



exp((t-t abs )^ cl )p(r), (2) 



with the classical Liouvillian S£ c \ for the full three-body 
Coulomb system propagated from the time i a b s of pho- 
toabsorption. The projector V ++ indicates that we have 
to integrate only over those parts of phase space that 
lead to double escape (the asymptotic final energies of 
the two electrons, e and E — e, are positive). The inte- 
gral in Eq. (^|) is evaluated with a standard Monte-Carlo 
technique which entails following classical trajectories in 
phase space. 

The electrons are described immediately after absorp- 
tion by the distribution 



(3) 



where JV is a normalization constant. The primary elec- 
tron absorbs the photon which has an energy fko. With 
S(fi) we demand the absorption to occur at the nucleus, 
an approximation which becomes exact in the limit of 
high photon energy JTT| . This approximation signifi- 
cantly reduces the initial phase space volume to be sam- 
pled. Regularized coordinates |l2|, [p| are used to avoid 
problems with electron trajectories starting at the nu- 
cleus (fi = 0). 

The function p 2 (r 2 ,p 2 ) describing the secondary elec- 
tron in Eq. (0) is given by 
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Pl(f 2 ,p 2 ) = W^,{r 2 ,f) 2 )b(£' 2 - £fl). 



(4) 
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FIG. 1: Photoionization double-to-single ratio. Circles: 
benchmark experimental data (Samson et al. [Q). Full line: 
complete theoretical result. Dashed line: knock-out mech- 
anism only. Chained line: shake-off mechanism only. The 
arrow indicates the asymptotic ratio (~ 1.645%). 



It is obtained by calculating the Wigner distribution, 
^>{t2,P2), of the orbital ^(ra) = *o(fi = 0, r 2 ) for a 
choice of initial wavefunction \&oj and restricting the ini- 
tial energy of the secondary electron, eif, to an energy 
shell £b- In the KO mechanism the initial state correla- 
tion is not important so we take the independent particle 
wavefunction ^ (ri,r 2 ) = (Z% s /w) exp(-Z cff (r\ + r 2 )) 
with effective charge Z e g = Z — 5/16. From this choice 
follows e B = -Zl s /2 and e 2 u = p 2 2 /2 - Z cS /r 2 . 

The double-to-single ratio in the absence of the SO 
mechanism is simply given by 



Rko = P K £/(l 



(5) 



In Fig. [l| we show Rko as a function of the excess en- 
ergy E (dashed line). The shape is characteristic of an 
impact ionization process |I3 . For high energies the pri- 
mary electron moves away so quickly that there is no time 
to transfer energy to the secondary electron, Rko thus 
drops to zero as expected. The non-zero asymptotic ratio 
(indicated by an arrow in Fig. Q) is due to SO which we 
describe next. 

In contrast to KO the shake mechanism is inherently 
non-classical in nature. Moreover, initial state correla- 
tions are important for shake-off. As a generalization of 
the standard formula for SO Aberg gave an expres- 
sion for the probability to find the shake electron in state 
(f> a at any excess energy H, 



with 



r(r 2 ) = / d 3 r 1 ^(r 1 )* (r 1 ,f 2 ), 



(6) 



(7) 



where v(r\) is the wavefunction of the primary electron 
after it has left the atom. If it was in an s-state be- 
fore the absorption it is in a p-state afterwards. The 
secondary (shake) electron does not change its angular 
momentum. It can be found with probability P a in an 



> 

0J 



T3 



a. 



0.9 
0.8 
0.7 
0.6 
0.9 
0.7 
0.5 
0.3 









6eV 


_ 















2 


4 


e 














21 eV 









FIG. 2: Differential probabilities for separate knock-out and 
shake-off mechanisms for a number of excess energies. Circles: 
knock-out mechanism results from binning. Solid lines: fits 
through circles. Dashed lines: shake-off mechanism results. 
See text for details. 



hydrogenic eigenstate of the bare nucleus, being cither 
bound (a — n 2 ), or in the continuum (a — e). As for KO 
we assume that the primary electron absorbs the photon 
at the nucleus. In this situation we do not need to know 
v{f{) but can simply replace ip v (r 2 ) by V?o(^i = 0, r^) 
in Eq. M). We may further simplify the calculation of 
shake-off for practical applications in two-electron atoms 
by taking for ^0(^1 = 0,r2) a normalized hydrogenic 
wavefunction 4>ff°(r 2 ) where the correlations have been 
'absorbed' into an effective shake charge Zso j^j- F° r 
Zso « 2 — 0.51 the exact asymptotic ratio Roo — 0.01645 
fig, [r^j is reproduced. We have found little difference for 
the shake probability as a function of excess energy be- 
tween this simple ansatz and a fully correlated Hylleraas 
wavefunction |Q for "Jo. 

The shake-off probability of Eq. (0) reduces now to 



|«>a|0 ZSO } 



(8) 



The total double ionization probability from shake-off at 
finite energies E is given by integrating expression (^) 
over the energy e of the shake electron in the continuum 

(a = e), 



p ++ 
^so 



(E) = f 
Jq 



deP, 



(9) 



The photoionization ratio when only the SO mecha- 
nism is taken into account (same as Eq. (||) but for SO) 
is shown in Fig. |Tj (chained line). The ratio rises slower 
than the KO mechanism result up to an energy of around 
100 eV where the KO ratio reaches its maximum value. 
The SO ratio continues to rise until at a couple of hun- 
dred eV it moves more slowly up toward the asymptotic 
value. An interesting feature of the plot is where the KO 
and SO results cross at an excess energy of ~ 350 eV. 
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FIG. 3: Singly differential cross sections normalized to 1 for 
e = 0. Solid lines: our complete theoretical results at excess 
energies of (a) 6 eV, (b) 21 eV, (c) 41 eV. Dashed lines: new 
results of Colgan et al. |^l[ at excess energies of (a) 4 eV, (b) 
20 eV and 25 eV, (c) 40 eV. 
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To obtain more insight into the two mechanisms we 
calculate the differential probabilities dP^" + /de, where 
X stands for either SO or KO. In our classical model of 
the KO mechanism we divide the interval of values for e 
which corresponds to double escape (0 < e < E) into N 
equally sized bins (we take N = 21) and work out the 
differential probability by finding the trajectories which 
fall into the bins. For the SO mechanism the probability 
per energy unit, P e in Eq. (^J), already gives the differen- 
tial probability. Since the electrons are indistinguishable 
the differential probabilities must be symmetrized about 
the equal energy sharing point (e = E — e = E/2), 



dP++ 



de 



dP++ 



(e,E) 



dP x +- 



(E — e, E) 



de 



de 



(10) 



In the case of low excess energy (6 eV) we find a slightly 
concave shape for the KO distribution, see Fig. ||. This 
implies a preference for equal energy sharing, the typical 
behavior close to threshold Jl9{ . The SO result in con- 
trast displays a slightly convex shape which becomes flat 
as E — > +0. Unequal energy sharing is always preferred 
by SO since the photoelectron is fast with respect to the 
secondary electron. For all higher excess energies shown 
both mechanisms display a convex form. 

SO may be viewed as an additional quantum contribu- 
tion to the quasi-classically calculated double photoion- 
ization given by KO. This means that the full result is 
given by 



da++ 
de 



— Cabs 



de 



aJ so 
de 



(11) 



Integration over e yields the total double ionization cross 
section, 



CT ++ = <7 abs (P+++P 



SO 



'KO t u SO • 



(12) 



FIG. 4: Absolute singly differential cross sections. Solid lines: 
our theoretical results at excess energies of (a) 6 eV, (b) 11 eV, 
(c) 21 eV, (d) 31 eV, (e) 41 eV. Circles: recalibrated (see text) 
experimental data of Wehlitz et al. |^2j at the same excess 
energies apart from (a) which is at 5 eV. The triangles in 
(e) additionally show the Wehlitz data renormalized to the 
cr++ (41 eV) of Samson @. 



The single ionization cross section is er + = <7 a b s — c ++ 
and the double-to-single ratio is given by R = <r ++ /a + = 
P++/(l - P++), where P++ = P++ + P++. 

In Fig. |l| we compare the ratio R (solid line) to the 
experimental data of Samson et al. [ p0| . For excess 
energies up to 200 eV we find an excellent agreement. 
In the energy regime where the two contributions are of 
the same size there is a deviation between experiment 
and our result (at worst 8%). Exactly in this situation 
any interference which exists between SO and KO would 
show its largest effect. The deviation we find may be 
due to such an interference which we cannot account for 
since we determine KO quasi-classically. At higher ener- 
gies the difference decreases again (already visible in the 
plot) until at very high energies our result reproduces the 
asymptotic ratio. 

Knowing that the differential probabilities for KO and 
SO enter the full ionization probability with the same 
weight we can assess the relative importance of both con- 
tributions at different energies. From Fig. |^ one sees that 
at 110 eV SO has become more important than KO for 
highly unequal energy sharings. As energy is increased 
to 450 eV SO begins to dominate regions of unequal en- 
ergy sharing. On the other hand, KO is higher at equal 
energy sharing for all excess energies E. 

Fig. H shows that our singly differential cross sections 
(SDCS) agree well with the recent ab initio theoretical 
results of Colgan et al. |2l|]. We note that the results 
of Proulx and Shakeshaft [E3| show a concave shape for 
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FIG. 5: P^o (solid line) as a function of the scaled energy 
E/E' B compared to the cross section for electron impact ion- 
ization of He + ^(| (circles) as a function of E/Eb (see text). 
Additionally the impact ionization data has been multiplied 
by a factor C = 4.67 x 10 15 cm -2 to make the maxima of 
both curves the same height. (1/C may be interpreted as the 
geometric cross section.) 



excess energies below 20 eV. This is in disagreement with 
our results which are convex down to 6 eV. In Fig. || we 
compare our absolute SDCS to the experimental data of 
Wehlitz et al. [^2| which has been recalibrated using the 
values of the photoabsorption cross section of Samson et 
al. In their original work Wehlitz et al. normalized 
their SDCS using a photoabsorption cross section ]24f| 
which is now known to overestimate the 5-41 eV range 
by 9 to 16%. In addition, the photoionization ratio they 
measured at 41 eV is ~ 17% higher than the Samson et 



al. j20j data and so we renormalize Fig. |4|(e) to take this 
into account. 

Our approach not only facilitates the calculation of 
double photoionization, it also offers considerable insight 
into the physical process, e.g., concerning the similarity 
with electron impact ionization of He + pq |. Indeed, we 
can show that impact ionization may be viewed as the 
KO part of double photoionization (Fig. |J). The only dif- 
ference is that impact ionization sees a He + hydrogenic 
target electron with binding energy E-q = —Z 2 /2,Z = 2, 
while the KO process involves a bound electron with en- 
ergy E' B — —Z^ s /2. One may thus say that both pro- 
cesses differ only slightly, namely in the energy scale set 
by the respective bound electron. 

We conclude that the separate formulation and cal- 
culation of knock-out and shake-off offers an accurate 
description of double photoionization. In principle this 
approach can be extended to describe angular differen- 
tial cross sections. We have used a description in terms 
of the simplest wavefunctions possible having in mind to 
tackle three electron problems in a similar way. 
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